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Abstract 

We propose a novel approach to solve the three-nucleon (3N) Faddeev equation which avoids the compli- 
cated singularity pattern going with the moving logarithmic singularities of the standard approach. In this 
new approach the treatment of the 3N Faddeev equation becomes essentially as simple as the treatment of 
the two-body Lippmann-Schwinger equation. Very good agreement of the new and old approaches in the 
application to nucleon-deuteron elastic scattering and the breakup reaction is found. 



PACS numbers: 21.45. +v, 24.70. +s, 25.10.+S, 25.40.Lw 



I. INTRODUCTION 

The three-body Faddeev equations in the continuum put three-body scattering on a firm 
mathematical basis. The complicated asymptotic boundary conditions in configuration space [3] 
are fully included. At the time of their formulation the computer power, however, was insufficient to 
solve them directly given local two-body forces. Therefore one started with forces of finite rank (the 
simplest being a rank one separable force), which turns the three-body problem in a partial wave 
basis into a finite set of coupled integral equations in one variable - and thus feasible at that time. 
Nevertheless the way it was formulated in momentum space the free three-body propagator lead 
to a complicated singularity structure of the integral kernel, namely to logarithmic singularities, 
whose position depend on the external variable - so called moving singularities [3]. In order to 
avoid this obstacle the path of integration in the kernel is moved into the complex plane (contour 
deformation), which imposes of course conditions on the analytical properties of the two-body 
separable forces. The need, however, for using the realistic two-body forces, which are dominantly 
local, enforced their expansion into a series of finite rank forces [4], which was tedious and finally 
overcome by integrating the logarithmic singularities directly on the real momentum axis [5] . It took 
until the eighties [6J that fully realistic two-body forces could be handled, now in a set of coupled 
integral equations in two variables. That approach was based on Spline interpolations Q], which 
allowed to integrate analytically over the interpolated Faddeev amplitudes under the integral. In 



reliable manner 



this manner in the three-nucleon system the realistic high precision nucleon-nucleon (NN) forces 
[9. [lo| together with three-nucleon forces (3NF) of various types [ll|, [ijj] can be handled in a fully 



These studies are a basic foundation for testing nuclear forces and triggered 



a tremendously rich set of neutron-deuteron (nd 
the world. A representative set can be found in 
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and proton-deuteron (pd) experiments all over 
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301 ] . Especially the need of 3N forces in conjunction with current two-nucleon (2N) 



forces was firmly established 
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32, 



33, 



34|. In recent years that interplay of 2N and 3N forces 



has been substantiated by the theoretical insights into nuclear forc es g ained through effective field 
theory for 7r's and nucleons (A's) constrained by chiral symmetry 35J. 

The feasibility of controlling the 3N continuum in the Faddeev scheme opened the door to 
evaluate reliably the final state (initial state) interaction in photon induced processes on 3 He (pd 
capture processes) 36J. This provides important insight into the electromagnetic nucleonic current 
operator, the interplay with 3N forces and into properties of electromagnetic nucleon form factors. 

Despite that achieved technical status of controlling the integral Faddeev kernel in the three- 
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body continuum it would be desirable to avoid those logarithmic singularities totally. A suggestion 
in that direction was undertaken in [3tJ where, however, the presence of the virtual state pole of 
the 2N t-matrix in the state 1 5o caused problems. It is the purpose of this paper to establish a 
definite solution of that long lasting technical challenge with the moving logarithmic singularities 
by formulating the three-body Faddeev kernel in the continuum in such a manner that only trivial 
poles occur in one variable, which can be handled as easily as in the 2-body Lippmann Schwinger 
equation. 

In section II we give a brief reminder of our standard approach followed by the discussion of all 
possible alternative choices of variables in the intermediate state integral of the Faddeev integral 
kernel. One of them sticks out which is free of all technical difficulties and will be displayed in 
detail. In section III we compare numerical results for 3N scattering in the old and new approach 
using modern high precision NN forces. The great simplification achieved with the novel approach 
shows up also in case of finite rank 2-body forces as displayed in section IV. We conclude in section 
V. 

II. DIFFERENT FORMS OF THE 3-BODY FADDEEV EQUATION IN THE CONTIN- 
UUM 

We use the notation detailed in The nucleon-deuteron (Nd) breakup amplitude is 

<pq\U Q \$ >=< pq\(l + P)T\$ > , (1) 

where p and q are standard Jacobi momenta, P the sum of a cyclical and anticyclical permutation 
of 3 particles and |<3? > the initial product state of a deuteron and a momentum eigenstate of the 
projectile nucleon. 

The amplitude T|<3? > obeys our standard Faddeev type equation 

T|$ >= tP\§ > +tPG T\<S> > , (2) 

where t is the 2N off-shell t-operator and Go the free 3N propagator. Introducing the momentum 
space 3N partial wave basis \pqa > and projecting ([2]) on these states we get 

(pqa\T\$) = (pqa\tP\<S>) + (pqa\tPG T\$) . (3) 

The kernel part (pqa\tPGoT\&) can be evaluated using the completeness of the states \pqa > 



MS 



{pqa\tPG T\$) = ^ J ' p l2 dp 'q ,2 dq ' p" 2 dp" q" 2 dq" (pqa\t\p 'q 'a ') 



Here we use the Balian-Berezin 



38, 



39] approach to calculate the permutation matrix element 



{p 1 q ' a'\P\p" q" a") (see Appendix for details) 



(pqa\P\ P 'qW) = /' dJ-^^ G™(q,q>, X ) . (5) 



_1 p 

Moreover the 2-body t-matrix conserves the spectator momentum q and all discrete quantum 
numbers except the orbital angular momentum I: 



<pqa\t\p'q'a'> = ^ (pp E(q) = E - —q 2 ) 

$s a s a ,5j a j a ,S tat A„A a ,Si a i a , • (6) 



Finally we extract the deuteron pole in the channels \a >= \a d > which contain the 2-body 
3 Si — 3 D\ states. Thus we define 

f s a j a t a( /. F( u _ K 3 C a (PiP'; E (q)) m 

t u& (p, P ,^))___^_ (7) 

for the deuteron quantum numbers s a = j a = l,t a = 0,l a ,la = 0,2 and keep t as it is otherwise. 
That pole property obviously carries over to the T-amplitude and we define just for the \ a >= \a d > 
channels 



(pqa\T\& 
ie 

Using all that the coupled set (J3J) is for a/aj 



(pw\Tm = - • (8) 

h + l€ - 4^1 - h d 



(pqa\T\$) = (pqa\tP\$) [ p ' 2(ip ' p" 2 d P" q" 2 d( l" 

l& a" J 

4 p ' p" 

(v"a"a"\f\<§>) - 1 

^ ^ e + ie - ^q" 2 - E d a a ^ 1 1 11 E + ie - ±{p" 2 + \q" 2 ) K) 

where 5 a " a " = 1 — d a "ot"i an( i f° r a = a d 

d d 

(pqa d \f\$) = {pqa d \tP\<5>) fp^P ' p" 2 d P q" 2 dq" 

la a" J 

1 dx trjf^p -E { q)) G r al ^ q \ X ) 5{p '7 l)5{p ''7 2) 

4 p p" 

{p"q"a"\f\<5>) 1 

d E + ie - -^q" 1 - E d d E + ie-± (p" 2 + ±q" 2 ) 
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The shifted p ' and p " arguments are 



n-j - \/ q" 2 + ^q 2 + qq"a 



= \jq 2 + ^q" 2 + qq"x . (11) 

We are left in Eqs. and (jlOp with four integrations over p ' ,p" , q" , and x, where any two of 
them can be performed analytically using the two <5-functions. Thus there are six possibilities, which 
we regard now in turn and we discuss the advantages or disadvantages to use them. Apparently 
it is sufficient to discuss only the kernel parts and moreover just one, say in ([U]), which we shall 
denote just by "the kernel". 

A. Analytical integration over p ' and p" 

The resulting form of "the kernel" is 

(pqa\tPG T\$) = EE / f dx t*£ ta (p,7r i; E{q)) 



BB , „ UJC (7r 2 q"a^\T\<S>) 

4m 

(12) 



d E + ie - £-q" 2 - E d d 



E + ie- ±(q 2 + q" 2 + qq"x) ' 

This is our standard approach [f], Q. For q < q max = yj %raE the integration over x leads to 
logarithmic singularities depending on q and q " - the so called moving singularities. Nevertheless 
the advantage is, that the complex g"-dependence of {iT2q"a"\T\<&) can be properly mapped out. 
That nontrivial dependence of the T-amplitude arises from the property of the 2N t-matrix in the 
state 1 S'o and from the 3N breakup threshold behavior 13]. 



The deuteron pole at q" = qo = y ^(E — Ed) can be taken care of in the {/'-integration using 
e.g. a subtraction method. 

B. Analytical integration over x and p " 

We rewrite the (^-functions in ([9]) and (fTUl) as follows 



with 



6(p'-n l ) = -^ F 8(x-x )@(l-\xo\) , (13) 

qq 



_ p' 2 - l/iq 2 - q" 2 _ p" 2 - l/4q" 2 - q 2 

X — 7. — 7. ) (.14 J 

qq" qq" 
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and 



5{p" - n 2 ) = S(p " - ^p '2 + If - fq «2 )e (p ' 2 + ^ 2 - \q " 2 ) . (15) 



4 

The resulting form of "the kernel" is 

(pqa\tPG T\S) = ~J2J2 [p' d P 'q ! 'dq''t s i :if-{p,p'-,E{ q ))G^, l {q, q '',x Q ) 

q la a" J 

l P ' 2 - l/4g 2 - g //2 |W /2 3 2 3 //2 
99" 4 9 ~4 9 



with 



P" = \/P' 2 + 1<1 2 ~ ~9 " 2 • (17) 

The two 0-functions define the domain D for the integrations over p ' and q" which is an open 
rectangular region in the p ' — q" plane restricted by the straight lines q" = | ± p ' and q" = p ' — | 
as displayed in Fig. [TJ 

The integrations over p ' and q" in (|16p split now in the following way 



rq/2+P ' (n"n"ry"\f\<S>\ 

/ JWG™„{qy,x )(6« la ,, - t [P Q J 2 ; +WWV / m$)), (18) 

J\q/z-p 'I jfc/ + ze_ 4^^ 

For the channels a" different from only a simple pole in the p' variable occurs, positioned at 

Po = \J\{qmax ~ 9 2 ) an d therefore of concern only for q < q m ax . At q = q m ax or po = there is no 

pole since the q" integral vanishes. For the channels a" = a'^ the q"- integral contains the deuteron 

pole at q" = qo = \J — E^). If q$ does not coincide with the limits of integration || — p '| and 

| -\-p 1 that integral generates a smooth function of p ' . In case it coincides, however, logarithmic 

singularities occur. Their positions are determined by 

,9 n 

90 = I , 

90 = \ +p' ■ (19) 

We can avoid them in the following manner. The product of the free propagator and the deuteron 
pole term can be written as 

1 1 



E + ie-Mp' 2 + \q 2 ) E + ie-£-q" 2 -E d 



{ E + ie-±{p> 2 + \q>) E + ie-^-E^Ed + y* ' (2 ° } 

where we used ()17p . This provides separation of the free propagator and the deuteron pole singu- 
larities. Thus for the a" = cx' d channels alone one has the contribution to "the kernel" 

±J2E / p'd P 'tt:^( P , P ';E( q )) / q"dq"G™(q,q",x ) 

1 1 </</Vj|T|$> 



y E + ie-^{p' 2 + \q 2 ) \E d \ + ±p»z 

' 1 <AVj|f|$>] (21) 



in 



E + ie- f^q" 2 - E d \E d \ + 

The first part has only the free propagator singularity in the p ' variable like the part for the 
channels a" different from a d in (|18p and in the second part we change the order of integrations. 
Thus (|2ip yields the contribution to "the kernel" from a" = ol' d channels 

1 



-EE/ p'd P 't s r jf«(p,p>;E( q )) 
9/2+P' .. „„ .. <rr,"a"rv" 



E + ie -Up ' 2 + \q> 



/ WGBfe/.V Pq dl ' 

J\a/2-t> 'I d 



1 „V/2 



1,/2-p '| d \E d \ + ^P 

oo i 



fl/ 2 +l" . , <r r)"a"n"\T\& ■> 

/ p 'dp >t s r jf«(p,p 'lE^G^qy ^^y^y • (22) 

J\q/2-q"\ d \E d \ + -p" 2 

Now in the second part the p '-integration yields a smooth function in q" and the deuteron 
singularity in the (/'-integral is a simple pole. 
Altogether "the kernel" flE]) is 

oo 2 



2 f°° 

< pqa\tPG T\<S> >= - E E / P > 



(p,p ';E(q)) . ., 



4 

rq/2+P ' <- n"a"n/"\f\<b ~> 

/ 9 d<7 G aa"{l^ ^oXV'a" <P q a |T|$ > +* Q // Q " 

./la/2-» 'I 



1 „//2 



'[g/2-p'| lAil + ^P 

?££ fV-" 



rg/2+q" 



rll^+l . ^ <n" rt" ni"\T\(b ^ 

/ , ^wjgg^g^.xo) ^ y^T'! • ( 23 ) 

■/[«/2-9"| |Ai| + 

It is the best form out of the six possible ones for the numerical performance. Also the complex 
q-dependence of the T-amplitude near q l! — qmax in the channels cv containg the ^Sq two-body 
component can be properly mapped out like in case III Al 
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C. Analytical integration over x and p ' 



We choose one <5-function in ^ as in (|13|) and the other one as 



*(P" " ^) = ^(P ' " \/p" 2 + j<?" 2 " j9 2 W 2 + ~<T " > 2 ) • (24) 



This leads to "the kernel" 



{pqa\tPG T\$) = "EE / ^ / it**"^ '' E (l)) G ia »(?,?", *o) 

q U a" J J 



- \P" 2 - \q" 2 - Q 2 \W 2 + ~/ 2 ~ ~/) 



(p"q"a"\T\<S>) - , „ 1 

d E + ie- ^q" 2 - E d d E + le - ±{p" 2 + f q" 2 ) 



with 



p ' = yV /2 + 3/4q" 2 - 3/4g 2 . (26) 

The two 0-functions define the domain D for the integrations over p" and q" which is an open 
rectangular region in the p" — q" plane restricted by straight lines q" = 2(q±p") and q" = 2p" — 2q 
and which leads to the integrals 



r, poo p2q+2q" 

(pqa\tPG T\t>) = -EE / d P"P" / dqW't s l t ata (p,p';E(q))G^4Q,Q"^o) 
q l?y Jo J\ 2 <i- 2 v"\ 

(v"a"a"\f\$) - 1 

(5 ^ q Vrf^ + 5 a ,, ol n(p" q "a"\Tm) f— ^— , (27) 

^ a d E + ie - ±q» 2 - E d a ^ 1 1 l! E + ie - ±(p" 2 + \q" 2 ) V ' 

One can proceed analogously like in III Bl and would arrive also at a form free of logarithmic 
singularities. The form III Bl might appear, however, more favorable since each propagator appears 
only with one integration variable. 



D. Analytical integration in q" and p ' 

We keep the first <5-function in ([9]) as it is and rewrite the second one as 

5(p" - vr 2 ) = - P -5(q" - (-2qx + 2^p" 2 - q 2 (l - x 2 ))9(p" 2 - q 2 ) . (28) 

yp 2 — q 2 [l — x 2 ) 

The resulting form of "the kernel" is 

tt± ta (p>^ E (v)) 



(pqa\tPG T\<5>) = f p"dp"®{p" 2 - q 2 ) f dx 



y / p" 2 -q 2 (l-x 2 ) 
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nBB i i, we (p"q"a"\T\<S>) - „„„ 

d £ -Me - ^g" 2 - ^ d 

' (29) 



£ + it - p" 2 - 3( v / p" 2 -g 2 (l-x 2 ) - gx) 2 
with 



g" = 2( V / p" 2 -q 2 (l -x 2 ) - qx) . (30) 

This form appears quite complicated in the pole structure of the free propagator as well as for 
the deuteron pole and the 1 Sq virtual state pole singularity close to q" = q ma x and we do not 
consider it further. 



E. Analytical integration over q" and p" 

We rewrite the first (^-function in {TOT) as 



p'S(q" - U p 12 - k 2 (l - x 2 ) - f )) 
S(p ' ~ TTi) = V 4 2 &(P ' 2 ~ 1/V) , (31) 



^p' 2 -\q 2 {l-x 2 ) 

and keep the second as 



S(p " - v / (q-l/2q") 2 + qq"{l + x)) . (32) 

Thus p " is always greater equal zero. It is zero for x = — 1 and p ' = 3/2q. 
"The kernel" is 



(pqa\tPG Q T\$) = ^^f* '6(p ' 2 - l/4g 2 

a" J 



dx t s «jf«(p,p >;E(q)) U G™„(q, q ",x) 



1 pl 2_l2 [l _ x 



,.2\ 



(v"a"u"\T\<&) - 1 

(<W „ ^ 9 3 ,2 F + ^WW)^ 1! ,2^3 2 , > (33) 

d E + ze - ^g" 2 - E d d E + it - ±(p ,z + fq 2 ) 



with 



and 



V " - \^-\q ! '? + qq"{i + x) , (34) 



P /2 -^ 2 (l-* 2 )-f • (35) 

Though the free propagator singularity is again a simple pole in p ' a two-fold integration in the 
T-amplitude is required. This appears to be a disadvantage (though surmountable) against the 
case III Bl However, in the deuteron pole both integration variables p' and x occur, which is quite 
unfavorable. 
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F. Analytical integration over x and q " 

The first ^-function in Q is taken as in (|13|) and the second one as 

S(p " - vV + l/4g" 2 + qq"x) = " - ^4/3(p > 2 - p" 2 ) + q 2 ) 

e(4/3(p' 2 -p" 2 ) + q 2 ) . (36) 

It results in "the kernel" 

(pqa\tPG T\$) = J dp ' p ' J ^V'e(4/3(p /2 - p" 2 ) + g 2 )9(l - |x |) 



la a 1 ' 
3a.t ai lt - l/iUrj BB 

J" (A Virjg) ,„„ „ |T ™ 

d E- E d + ie- ±(p' -p" 2 )-^q 2 d 

' (37) 



E + ie-Up' 2 + \q 2 ) ' 



with 



q" = \Jl(p /2 -p" 2 ) + q 2 ■ (38) 

Though the free propagator singularity is just a simple pole in one variable in the deuteron pole 
both integration variables p ' and p " occur, which is less favorable. 

We conclude that case MB I is clearly the most favorable choice and we compare in the next 
section results for different 3N observables choosing our standard approach, case III Al and that 
new one, case III Bl 



III. COMPARISON OF THE NEW (|H"Bj) AND STANDARD (UTAl) APPROACHES 

Since the complicated singularity pattern in the old approach exists only for q" < q ma x we 
applied the new approach only there and kept the old one for q" > q ma x, where only a simple 
deuteron pole is present. We used the CD Bonn [9] potential restricted to act in the two-nucleon 
partial wave states with total angular momentum j < 1. In Fig. [2] we show the resulting nd elastic 
scattering angular distributions for the cross section and various analyzing powers at an incoming 
neutron lab. energy El l . = 13 MeV. The agreement obtained with the two approaches is very 
good. The cross sections at the same energy for two geometries of the Nd breakup are shown in 
Fig. [3l Again the agreement is very good. 
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IV. FINITE RANK FORCES 



Using the choice ITTBl simplifies the treatment of the 3N continuum in case of finite rank 2-body 
forces also very significantly since there are no longer logarithmic singularities. For the sake of a 
simple notation we keep only s- waves and restrict the 2-body force to act only in the states 1 5o 
and 3 Si. This leads to two coupled equations for the two amplitudes Tk(pq) where k = 1 (k = 2) 
goes with 1 Sq ( 3 #i), respectively. Choosing the kernel of the type ([23]) one obtains explicitely 



2 f°° 

Txipq) = T?(pq) + - / dp 'p 'ijfop '; E{q))- 
Q Jo 



1 



E + ie-±(p' 2 + h 2 ) 
q/ 2 +p ' fo(n"a") 

.9/2-* 'I E d + ^P" 2 

2 r°° 1 ri/ 2 +i" fi(v"a"') 

qJo E + te - ^q" 2 - E d J\ q /2- q "\ E d + ±p" 2 

2 1 

T 2 (pq) = T°(pq) + - / dp 'p %(p,p '; E(q))-—, 

QJo E + ie- ±(p 12 + %q 2 ) 

[il 2+ P ' fo(v"a") 
/ cW'(G 21 WY) + G 22 ^ff4) 

■/|9/2-p '| #d + „P" 2 

2 f°° 1 /•9/2+?" fi(n"a") 

~ ~ / <W'-— 3 „ 2 F g 22 / dp ' P 't 2 (p,P ^(g)) -^ "I (40) 

g Jo £ + ie - ^q" 2 - E d J\q/2- q "\ E d + ±p" 2 

with 



//'. x / p /2 + 3 g2 _ 3 ?//2 (41) 



Now we assume the finite rank forms 



h{pp';E(q)) = gi(jp)ri{E(q))gi(jp ') , 

* 2 (PP ';£(?)) = g 2 (p) F ^. f2 ™l t-i 92(P ') , (42) 

e + ^ - i^q -Ed 



and obtain 



2 f 00 1 

Txipq) = T°(pq) + gi{p)n(E{q))[- / dp 'p ' 9l (p ')— 

9 Jo £ + *e-^(p' 2 + fg 2 ) 

/•9/2+p ' fi(v"a") 
/ WiGuTtfft + Gu j; 2 ** /2 ) 

•Zk/2-p 'l + ib" 2 

2 f°° 1 rq/ 2 +q" fo(n"a") 

~ l dq"q"— \ - - G 12 dp 'p ' gi {p ') * 2(P ? \ ] (43) 



qJo E + ie-^q" 2 -E d ^ J\ q /2- q »\ ' " * ^ E d + ±p" 2 

2 f°° 1 

T 2 (pg) = T$(pq) + g 2 (p)f 2 (E(q))[- / dp 'p ' g 2 (p ')——, 

9 Jo E + ie - ±(p ' 2 + ±q 2 

" q/2+p ' f 2 {p"q") 



rq/z+p 

/ dq"q"(G 21 T l (p"q")+G 22 
J\ q /2-p '| + ±P" 2 



Q Jo 



oo 1 rq/2+q" f t I, < ll\ 

Wt— 3 ,n r G22 / dp ' p ' 92{p H I I J (44) 

E + ie- j^q" 2 - E d J\ q /2- q "\ E d + -p" 2 
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where going through the same steps it results 



jf(pq) = Ng^niEiq^iq) , 

f°(pq) = Ng 2 (p)T 2 {E{q))F${q) . (45) 

The normalisation factor N provides the dependence on spin and isospin quantum numbers (m d , 
mo are spin projections of the initial deuteron and nucleon, respectively, and the nucleon's 
isospin projection), on the initial momentum q$ and the deuteron normalisation factor N d defined 
as UP) = N d f%: 

Efj 

N = 4=*M TW (l^,m(i,mo,M)— . (46) 
V47T 2 1 q 



Further F?(q) and F^(q) are given as 



2 /-90+9/2 , , , , ,M\l p 12 + h 2 - ho) 

12 i 3„2_3„2 



9 /*<20+<?/ 2 ( 

= Gi 2 - / dp'p' gi {p')- 

9 7|<7n-<7/2| 



9 y| m -,/2| Ed - P ^ 

a m 

= G 22 - / dp >p 'g 2 (p ') V * 3 2 



It follows the structures 



(47) 



Tiipq) = giWniEiqfiF^q) , 

f 2 (pq) = g 2 (p)T 2 (E(q))F 2 (q) , (48) 

and therefore one obtains the two coupled one-dimensional equations inserting explicitely the in- 
tegration limits 



2 f°° 1 

r' 2+r ' m lM E«-m«>) + 

qJo E + ie- ^q" 2 - E d J\ q/2 - q »\ E d + ±p"< 

2 1 

^k/2-p 'I ^ + ^P" 2 

2 r , „ „ g 22 r /2+q " rl i i f r Mf^w))' 
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It remains to provide the factors 

G\i = -3- = G22 , 

o 

G12 = ~^ = G 2 i . (51) 

Note for E = -^q 2 there is no singularity at p ' = since the (/"-integral vanishes at p ' = 0. 
Also there occur only simple poles, which can be treated by subtraction. Using for instance Spline 
interpolation for Fi(q") and -^2 based on a set of grid points the two integrals can be trivially 
performed and one obtains a low dimensional inhomogeneous algebraic set of equations. 



V. SUMMARY AND CONCLUSIONS 



Starting from a partial wave decomposed form of the 3N Faddeev equation for a breakup 
amplitude we discussed the six possible choices of integrating over internal angular and momentum 
variables in the integral kernel. While our standard approach integrates over the moving logarithmic 
singularities along the real momentum axis we found a new one, which totally avoids that technical 
obstacle. It is very simple. The free 3N propagator singularity appears as a pole in a single variable, 
which can be taken care of trivially (like in the 2-body Lippmann-Schwinger equation). The 
deuteron pole singularity is again a simple pole and moreover the two poles are cleanly separated 
in two different integration variables. 

The other four choices for internal integration variables turned out to be less favorable. 

We numerically compared the two approaches evaluating some nd elastic and breakup observ- 
ables and found very good agreement. This should open the door to handle the 3N continuum as 
simply as solving the 2-body Lippmann-Schwinger equation - though of course some more variables 
occur. It will also simplify the application to electromagnetic processes in the 3N system, where 
initial and final state interactions have to be treated properly. 

We also draw attention to the Balian-Berezin method, which was proposed long time ago, and 
which deserves much more attention than it received up to now. 

The inclusion of 3N forces do not change the singularity structure of the kernel and can be 
equally well treated in the new approach but is left to a forthcomig study. 
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APPENDIX A: PERMUTATION OPERATOR 



In view of using the Balian-Berezin method 38|, l39(] for evaluation of the permutation matrix 
element {p ' q ' a'\P\p" q" a") it is adequate to change from jj coupling to LS coupling 



\ pqa >= ^ \P<lP >< P\<* > > 

P 



where 



< (}\ a >=< (lX)L(s-)S(LS)J(t-)T\(l S )j(X-)I(jI)J(t-)T >= yJjILS { 



The permutation matrix element 



is 13 



(Al) 



I s j 

X \ I ) ■ (A2) 
L S J 



(p >q '(3 '\P\p "q "f3 ") = 25 s 'S »<5t >t "h >L »S, 



\ S' s 



(-)* "ViT" 



1 i t ' 

2 2 

1 



T t 



<p'q '(I 'X ')L 'fi\P 12 P 23 \p "q "(I "X ")L '/X > 



Since the momentum space matrix element in (|A3|) is independent of \i one can put 



(p 'q'(l 'X ')L'\P 12 P 23 \p"q"(l "X")L ') 

^ J> 'Q '(l ^ ')L VlftiflBlP \ "(I "A ")L V> 



2L 



Further we use 



,, .,, _ 6(p'-p)6(q'-q) Lfl , , 

2 2 /a \p q ) i 



< p'q'\pq(lX)L[x >- 



p* q* 

and the linear relations among the Jacobi momenta of different types and obtain 



(p'q'(l 'X >)L>\P 12 P 23 \p "q"{l "X ")L') 

rn;£/4p > (P '« ')Y^n (P "<7 ") 



2L 



<J(p '-1/2$'- q ")5(p " + q'+ l/2q ") . 



(A3) 



(A4) 



(A5) 



(A6) 

Here enters the idea of Balian-Berezin [38[ • The sum over the products of spherical harmonics 
is a scalar and depends only on the scalar products among the four unit vectors. Due to the two 
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5- functions all those scalar products are fixed by the magnitudes of the four momenta p', q',p", and 
q" . Thus we introduce 



X(P 'q 'P "q ") = £ (p 'q (p "q ") 



(A7) 



Consequently X can be taken out of the integral and one obtains 



(p'q'(l'X')L'\P 12 P 23 \p"q"(l"X")L') = 



8vr 



/2 



n //2 



(A8) 



p p" 
where x = q ' • q. 

Finally we combine geometrical factors and write the permutation matrix element in the form 



(pqa\ P \p'q'a') = J dx 



p* 



P 



(A9) 



with 



tti = \l q' 2 + ^<? 2 + qq'% 



"2 = \/q 2 + ^<?' 2 + gg'z , 



(A10) 



and 



G^to'x) = (4vr) 3 / 2 s TT ,s MTMT ,^jhfx^j >i >§ >i>(-iy '+* ' i 1/2 1/2 * 

1/2 T £ ' 



LS 



I 


s 


il 




f 

I ' 


s' 


'1 




< A 


1 

2 




.. 


X' 


l 

2 


/' 


lj 


L 

V 


5 


J 

) 




L 


s 


J 

) 





1/2 1/2 s , v -, , , 

> >^ (l\L,miOmi)(l X L,mi'm\>mi) 

1/2 5 s ' 



(-r^- roj (p)^' ro .,(p')n'mw(? ') • 



(All) 



We use our standard notation / = 2/ + 1. It is assumed that the z-axis is along q and the 
momentum q' lies in the x-z plane. That leads to the following components of the q, q', p, and 
p ' vectors 



q' 



P 



p> 



[0,0, q] , 

[qWl-x 2 ,0,q'x] , 
[q'yjl-x^^q'x + hq] , 
[-l q 'y/l-x 2 ,0,-q-l 



q'x] 



(A12) 
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FIG. 2: The nd elastic scattering angular distribution and analyzing powers at an incoming neutron lab. 
energy Ef ah — 13 MeV. The solid line is the CD Bonn potential prediction using our standard approach 
of handling the logarithmic singularities. The dotted line is obtained with the new approach without 
logarithmic singularities. The 2-nucleon states are kept up to j m ax = 1 
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FIG. 3: Cross sections for the exclusive d(n,nn)p breakup at an incoming neutron lab. energy E™ ah = 
13 MeV. Curves as in Fig 2. The upper configuration is the final-state-interaction (FSI) geometry with 
the polar angles of the detected neutrons 9\ = 39°, 02 = 62.5° and the azimuthal angle tfiu — 180°. The 
lower configuration is the quasi-free-scattering (QFS) geometry with the polar angles of detected neutrons 
8i = $2 = 39° and the azimuthal angle (f>±2 — 180°. 
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